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13.  ASS  TRACT  (Muatnutn  200  tworOkJ 

Surface  second  harmonic  generation  (SSHG)  has  been  used  to  study  the  chemical 
and  physical  properties  of  liquid  interfaces.  Information  about  the  thermodynamic 
of  neat  liquid  surfaces  has  been  found  by  observing  the  SH  signal  as  a  function  of 
temperature .  Extensive  studies  using  monolayers  at  the  air/water  interface  have 
been  undertaken.  Studies  of  acid/base  equilibrium  and  monolayer  phase  transitions 
have  been  performed.  The  gas-liquid  coexistence  region  and  an  orientational  phase 
transition  has  been  observed.  SHG  was  used  to  .*■  tudy  the  interface  properties  of 
fused  silica/electrolyte  solution.  Orientational  relaxation  has  been  observed  on 
picosecond  time  scale  at  the  air/water  interface  by  time-resolved  second  harmonic 
generation.  The  data  suggests  that  the  ground  ar.d  excited  states  of  Rhodam.ine  6G 
have  different  equilibrium  orientational  distributions. 

Femtosecond  studies  of  electrons  ir.  liquids  have  characterized  the  geminate 
recombination  and  solvation  dynamics  of  electrons  in  polar  and  nonpolar  liquids. 
Photodetuchment  of  electrons  from  halide  ions  in  aqueous  solution  has  also  been 
studies.  Picosecond  studies  of  chemical  intermediates  have  examined  the  role  of 
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Abstract 


Surface  second  harmonic  generation  (SSHG)  has  been  used  to  study 
the  chemical  and  physical  properties  of  liquid  interfaces.  Information  about  the 
thermodynamics  of  neat  liquid  surfaces  has  been  found  by  observing  the  SH 
signal  as  a  function  of  temperature.  Extensive  studies  using  monolayers  at  the 
air/water  interface  have  been  undertaken.  Studies  of  acid/base  equilibrium 
and  monolayer  phase  transitions  have  been  performed.  The  gas-liquid 
coexistance  region  and  an  orientational  phase  transition  has  been  observed. 
SHG  was  used  to  study  the  interface  properties  of  fused  silica/electrolyte 
solution.  Orientational  relaxation  has  been  observed  on  the  picosecond  time 
scale  at  the  air/water  interface  by  time-resolved  second  harmonic  generation. 
The  data  suggests  that  the  ground  and  excited  states  of  Rhodamine  6G  have 
different  equilibrium  orientational  distributions. 

Femtosecond  studies  of  electrons  in  liquids  have  characterized  the 
geminate  recombination  and  solvation  dynamics  of  electrons  in  polar  and  non¬ 
polar  liquids.  Phctodetachment  of  electrons  from  halide  ions  in  aqueous 
solution  has  also  been  studied.  Picosecond  studies  of  chemical  intermediates 
have  examined  the  role  of  restricted  geometries  on  elementary  chemical 
processes. 
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1.  Introduction  to  Surface  Second  Harmonic  Generation 

1.1  theory  Optical  nonlinearities  occur  in  media  exposed  to  very  high  light 
fluxes,  such  as  those  associated  with  lasers.  One  of  the  nonlinear  optical 
effects  is  second  harmonic  generation  (SHG).  The  physical  interpretation  of 
SHG  is  that  a  laser  beam  at  frequency  co  interacts  with  a  nonlinear  medium  and 
generates  a  nonlinear  polarization.  The  latter  being  a  collection  of  oscillating 
dipoles  acts  as  a  source  of  radiation  at  2co  [1  ].  The  key  feature  of  SHG  that 
makes  this  experiment  feasible  is  that  SHG  is  electric  dipole  forbidden  in  the 
bulk  but  allowed  at  the  interface  where  the  symmetry  is  broken  [2,3].  Thus  SHG 
can  be  used  to  probe  selectively  the  interface  layer. 

The  second  harmonic  polarization  is  proportional  to  the  second-order 
nonlinear  susceptibility  of  the  surface  XS(2), 


Pi  (2<o)  =Xsyk(2>  Ej(co)  Ek(co)  (eq  1 .1 ) 

In  the  experiment,  the  specific  elements  of  the  surface  nonlinear  susceptibility 
can  be  measured,  which  are  directly  related  to  the  population  of  the  molecules 
at  the  surface,  NSl  and  the  nonlinear  polarizability  of  the  molecule,  a(2) 
averaging  over  the  orientational  distribution  at  the  surface. 

Xjjk  =  Ns  <  aijk(2)  >  (eq  1.2) 

Therefore,  second  harmonic  generation  can  give  information  about  molecular 
population  and  molecular  orientation  at  the  surface. 

1.2  experimental  method  The  experimental  arrangement  that  we  use  in 
our  SH  studies  is  schematically  shown  below. 
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Figure  1 

Surface  Second  Harmonic  Experimental  Setup 


The  magnitude  of  the  second  harmonic  signal  is  proportional  to  the  square  of 
the  population  of  species  at  the  interface.  By  measuring  the  polarization  of  the 
second  harmonic  light,  the  orientation  of  the  molecules  at  the  interface  can  be 
infurred. 


2.  Interfacial  Studies  of  Neat  Liquids 

2.1  The  Role  of  Hydrogen  Bonding  on  the  Orientation  of  Molecules 
at  the  Neat  Air/Liquid  Interface  Expanding  on  our  results  obtained  for 
neat  water,  we  examined  other  neat  solvent  systems  in  order  to  sort  out  the 
relative  contribution  to  the  orienting  force  made  by  the  dipole-dipole  interaction, 
ionic  interaction,  hydrogen  bonding  effects  as  well  as  entropy  effects  due  to  the 
presence  of  a  hydrocarbon  chain.  Assuming  an  "up"  versus  "down"  model,  the 
nonlinear  susceptibility  could  be  related  to  the  net  orientation  of  the 
molecules.[4] 

In  examining  a  series  of  neat  liquids  we  found  that,  in  terms  of  the 
temperature  effects,  that  the  solvents  could  be  grouped  into  two  broad  classes 
(hydrogen-bonding  and  non-hydrogen  bonding  solvents).  We  found  for 
hydrogen  bonding  liquids  that  the  SH  signal  exhibited  a  strong  temperature 
dependence,  where  the  surface  sensitive  signal  decreased  with  increasing 
temperature.  This  is  to  be  contrasted  with  non-hydrogen  bonding  liquids,  where 
there  was  a  negligible  temperature  effect.  Our  results  are  easily  explained  by 
considering  that  there  is  a  strong  orientating  field  of  the  surface  molecules,  due 
to  the  effects  of  hydrogen  bonding.  When  temperature  rises,  more  and  more 
hydrogen  bonds  break,  diminishing  the  orienting  force,  which  thus  causes  the 
SH  signal  intensity  to  drop. 

Hydrophobic  effects  was  also  found  to  strongly  affect  molecular 
orientation  at  the  air/liquid  interface.  Specifically,  our  SHG  studies  on  a 


AFOSR-88-0014 


homologous  series  of  alcohols  show  that  the  hydrocarbon  chain  has  a  large 
influence  on  the  interfacial  molecular  orientation,  as  evidenced  by  the  SH 
signal  enhancement  obtained  when  switching  between  water  and  methanol, 
while  from  methanol  to  octanol  the  SH  signal  does  not  change  significantly.  If 
one  assumes  a  simple  "up"  versus  "down"  orientation  model,  we  found  that  only 
one  methylene,  i.e.,  methanol,  is  enough  to  favor  one  orientation  to  the  other. 
Namely,  that  the  hydrocarbon  portion  of  the  molecule  tends  to  favor  the  "up" 
orientation  (toward  the  air).  This  important  result  is  also  independently 
confirmed  by  surface  tension  measurements. 

2.2  Energetics  of  Orientation  of  Water  Molecules  at  the  Vapor/Neat 
Water  Interface  In  an  earlier  study  we  concluded  that  the  preferred 
orientation  of  water  molecules  at  the  interface  is  the  protons  toward  the  bulk 
water.  To  study  the  energetics  of  the  preferred  orientation  of  "up"  vs.  "down" 
we  measured  the  temperature  dependence  of  the  SH  signal.  Using  a  simple 
model  of  the  orienting  forces  at  the  interface,  similar  to  that  of  a  system  of 
pernament  dipoles  in  an  external  field,  we  obtained  an  energy  of  orientation  of 
280  ±  80  cal/mole  at  295  °K,  i.  e.  about  1/2  kT  at  room  temperature. 

3.  Studies  of  Monolayers  at  the  Air/Water  Interface  Using  Second 

Harmonic  Gsr.crstion 

3.1  Effects  of  Long  Chain  Molecules  on  the  Population  and  the 
Orientation  of  Organic  Molecules  at  the  Air/Water  Interface 

Determining  the  nature  and  the  effect  of  cc-adsorptior*  of  chornicoi'y  Hi'forprit 
solute  molecules  at  the  air/liquid  interface  is  a  fundamental  issue.  It  also  has 
direct  and  immediate  impact  on  our  ability  to  understand  how  the  presence  of 
one  adsorbed  molecule  can  influence  the  interfacial  population  and  orientation 
of  another.  In  order  to  address  this  issue,  we  have  considered  how  a  lipid 
monolayer  can  influence  the  ability  of  a  molecule  in  the  aqueous  subphase  to 
"penetrate"  into  the  lipid-rich  liquid/air  interface  region. [5]  We  were  further 
motivated  to  consider  this  particular  study  since  it  has  direct  bearing  on  the 
nature  of  transport  across  membranes  and  bilayers,  as  well  as  on  our  interest  in 
developing  novel  two-dimensional  materials.  Previous  studies  have  been 
hampered  in  yielding  detailed  information  on  this  problem  since  the  techniques 
available  (such  as  surface  tension  measurements)  could  not  discriminate 
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between  the  different  solute  molecules  in  the  interface.  With  the  SHG  method 
we  do  have  the  resolving  power  to  discriminate  and  determine  what  is  at  the 
interface. 

It  was  found  that  a  fatty  acid  monolayer  has  a  profound  influence  on  the 
penetrating  ability  of  a  solute.  Using  p-nitrophenol  (PNP)  as  the  surface-active 
solute  we  saw  that  its  surface  population  is  dramatically  enhanced  as  the  lipid 
density  is  increased  (reaching  a  maximum  at  approximately  50  A2  per 
molecule  for  the  lipid  density).  This  was  found  when  either  pentadecanoic  acid 
(a  15  carbon  chain  acid)  or  hexadecanoic  acid  (a  16  carbon  chain  acid)  was 
used  as  the  insoluble  fatty  acid  lipid.  We  interpret  this  enhancement  in  the 


OH 


no2 


Figure  3.1 

p-NItrophenol  (PNP) 

PNP  population  as  resulting  from  the  preferential  solubilization  of  PNP  by  the 
lipid  in  the  lipid  rich  interface.  It  is  pointed  out  that  at  still  higher  lipid  densities 
we  see  that  there  is  a  reversal  in  the  PNP  population  dependence  on  lipid 
surface  density.  The  turn-over  occurs  when  the  fatty  acid  lipid  density  exceeds 
about  50  A2/molecule.  From  independent  measurements  we  know  that  the  fatty 
ac'd  monolayer  becomes  close  packed  (more  ordered)  when  its  surface  density 
is  increased.  The  critical  area  of  50  A2/molecule  for  the  lipid  above  which  the 
PNP  molecule  population  decreases  is  due  to  the  fact  that  the  available  surface 
area  is  smaller  than  is  needed  for  which  a  PNP  molecule  can  reside.  In  other 
words,  as  the  lipid  monolayer  becomes  close  packed  it  "squeezes"  out  the 
moderately  sized  PNP  molecule  from  the  surface.  In  light  of  the  dramatic 
influence  of  the  lipid  monolayer  on  the  composition  of  the  co-adsorbed  solute 
PNP,  we  expect  that  the  monolayer  should  also  influence  the  orientation  of  the 
solute  as  well.  In  confirmation  of  this  idea,  we  found  that  at  low  monolayer 
densities  (in  the  limit  for  which  there  is  an  enhancement  in  co-solute  population) 
the  PNP  orientation  does  not  change  significantly  from  its  nominal  value  of  60- 
63°.  Above  the  critical  lipid  density  of  50  A2/molecule  (in  the  limit  where  the  lipid 
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is  close  packed  and  there  is  concomitant  diminution  in  the  PNP  population)  we 
fin j  that  the  molecular  orientation  of  PNP  flips  up  more  toward  the  surface 
normal.  At  25  A2  per  molecule  we  determined  the  molecular  angle  of  PNP  to 
be  51°. 

3.2  Measurement  of  pH  and  pKa  at  the  Air/Water  Interface  We  want 
to  investigate  the  effect  of  asymmetric  forces  at  the  interface  on  an  equilibrium 
constant  like  pKa.  SHG  is  used  to  measure  the  relative  population  of  the  acid 

and  base  at  the  interface.  Consider  the  acid-base  equilibrium  at  the  surface  for 
the  insoluble  material  p-hexadecyl  aniline, 


C16H33 

Ka 

+  H20  ~ - 

nh3+ 

AH+ 

Figure  3.2 

Acid-base  equilibrium  of  p-hexadecyl  aniline 

s 

The  pKa  for  the  above  equilibrium, 

pKa  =  PH  +  log— na-  (eq  3,1  J 


C16H33 

+  h30+ 

nh2 

A 


s  denotes  the  properties  at  the  surface  and  Na  and  Nha+  are  the  number 
densities  (molecules/cm^)  of  aniline  and  anilinium  at  the  surface  respectively 
To  obtain  the  pKa  one  needs  to  know  apart  from  the  concentration  of  the  acid 

and  base  forms,  the  hydronium  ion  concentration  at  the  surface. The  surface  pH 
is  related  to  bulk  pH  through  the  surface  potential  O  by  the  following  relation, 

s  b  e<J>  _ 

pH  =pH  +  2  3RT  3.2) 


The  surface  potential  O  can  be  determined  using  the  Gouy  Chapman  model  of 
electrical  double  layers  [6,7]  as 


AFOSR-88-OOI 4 


O  =  sinh'^  eNHA+ 

V 


2CekT 


(eq  3.3) 


where  C  is  the  total  number  of  ions/cm3  in  bulk,  e  the  dielectric  constant  of  the 
medium,  T  the  temperature  in  kelvin,  k  the  Boltzmann  constant  and  e  the 
electronic  charge.  To  calculate  the  surface  potential  from  the  Guoy-Chapman 
model  the  surface  population  of  the  charged  form  (acid  or  base),  obtained  from 
the  SH  measurement  is  used.  We  thus  obtain  for  such  a  system  a  total  SH 
electric  field 

_  nA  _  iA  NHA+  _  ,  Na  _  jA  ,  ^  na  _ 

Etot=  ^ EAe'f  +  EHA+ =  f^EAe'<t>  +  (1  )  EHA+ 

(eq  3.4) 


where  Na+Nha+=  No  is  the  total  number  of  molecules  at  the  surface  which  is 
known  by  virtue  of  spreading  known  concentration  solution.  Ea  is  the  SH 
electric  field  of  aniline  at  a  sufficiently  high  pH  where  Na  =  No,  similarly  Eha+  is 
that  of  anilinium  at  a  sufficiently  low  pH  where  Nha+  =  No,  and  d  is  the  phase 
difference  between  SH  electric  field  of  the  two  species. 

Note  that  the  measurement  of  the  SH  signal  Etot  vs.  bulk  pH  can  be 
Na  Nha+ 

transformed  into  ^  or  vs.  bulk  pH,  since  Etot,  Ea,  Eha+  and  b  are 

measured  independently  using  SHG.  The  data  are  fitted  to  eqs.  1,2, 3, 4 

s 

allowing  the  pKa  to  be  the  only  adjustable  parameter.  The  resulting  fit  is  shown 

as  a  solid  line  in  figure  3.3.  The  value  of  pKa  obtained  in  this  way  is  3.6±0.2. 
pKa  in  the  bulk  for  a  long  chain  aniline  is  5.3  units  at  ionic  strength  of  1 N.  Thus 
a  difference  of  1.7units  in  pKa  corresponding  to  fifty  times  increase  in  the  acidity 
is  observed.  This  increase  in  acidity  arises  from  the  higher  free  energy  of  a 
charged  species  at  the  air/water  interface,  shifting  the  anilinium  acid-base 
equilibrium  towards  the  energetically  more  favorable  neutral  form.  This  is  a 
direct  evidence  that  the  chemical  properties  such  as  equilibrium  constants  and 
pH  at  surface  can  differ  significantly. 
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pH  in  the  bulk 

Figure  3.3 

SH  electric  field  as  a  function  of  bulk  pH  at  constant  ionic  strength  of  1M  .  Solid 
curve  represents  the  fit  using  Gouy-Chapman  model. 


3.3  Study  of  surface  diffusion  by  SH  fluctuations  Insoluble  film 
(monolayer)  at  the  air/water  interface  is  a  typical  two  dimensional  system. 
Analogous  to  three  dimensional  systems,  these  insoluble  films  exist  in  solid, 
liquid  and  gas  phases  depending  on  the  density  of  the  molecules  at  the  surface. 
However,  the  gas-liquid  coexistence  region  is  of  particular  interest  with  regard 
to  the  nature  of  phase  transition,  structure  of  the  coexisting  phases  and  the 
dynamics  of  the  motion  at  the  surface  of  the  different  phases.  This  coexistence 
region  has  been  studied  extensively  using  a  wide  variety  of  techniques,  but  the 
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results  from  different  experiments  are  not  in  agreement  with  each  other  resulting 
in  a  lot  of  controversies.[8-1 1]  Hence  we  want  to  use  SHG  to  investigate 
surface  heterogeniety  in  the  gas-liquid  coexistence  region. 

A  heterogeneous  surface  in  the  gas-liquid  region  consists  of  two  different 
types  of  molecular  arrangements.  In  one  configuration,  the  molecules  are 
clustered  together  at  a  density  corresponding  to  a  liquid-like  surface  phase.  In 
the  other  arrangement,  the  molecules  are  sufficiently  far  apart  so  that  there  are 
no  interactions  between  individual  molecules;  this  corresponds  to  a  gas-like 
phase  of  the  molecules  at  the  surface.  In  the  coexistence  region,  the  surface  is 
heterogeneous  in  having  regions  of  liquid  clusters  and  regions  of  gas, 
interspersed  among  each  other.  The  liquid  phase  gives  rise  to  high  signal 
whereas  gas  phase  gives  small  signal  and  hence  the  time  dependent 
fluctuations  in  the  SH  signal  is  used  to  investigate  the  motion  of  these  liquid 
clusters  into  and  out  of  the  probe  area.  The  system  we  have  studied  is 
hexadecyl  aniline  spread  on  water  at  coverages  corresponding  to  the  gas-liquid 
coexistence  region.  The  surface  pressure  vs  area  phase  diagram  is  shown  in 
figure  3.4,  next  page. 

The  measured  SH  signal  l(t)  at  time  t  is  related  to  A|  and  Ag,  the  fraction 
of  the  area  under  observation  covered  by  liquid  and  gas  phases  respectively, 
as 


l(t)  =  l|A|  +  IgAg  =  ( 1 1 - lg )  A|  +  lg  (eq  3.5) 

l|  and  lg  are  the  signal  from  homogeneous  liquid  and  gas  phase  respectively. 
We  then  define  Al(t)  =l(t)  -  <!>,  where  <l>  is  the  average  signal,  and  obtain  a 
correlation  function, 

<AI(0)  Al(t)>  =  (l|-lg)2  <AA|(0)  AA|(t)>  (eq  3.6) 
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Figure  3.4 

Phase  Diagram  for  4  monolayer  systems 
1  Cl  6-phenol/water;  2  C-16-anallne/water; 

3  C-18-phenolate/2N  KOH;  4  Cl6-anlllnium/iN  HCI 


The  time  dependence  arises  from  the  changes  in  the  fraction  of  the  area 
covered  by  the  liquid  phase.  Assuming  the  cluster  size  to  be  invariant,  the 
changes  can  be  attributed  to  the  number  density  fluctuations  of  the  liquid 
clusters. 

The  characteristic  time  of  decay  of  the  correlation  function  is  thus  related 
to  the  time  that  a  cluster  of  liquid  remains  within  the  beam.  If  the  cluster  motion 
is  diffusive,  a  diffusion  constant  may  be  estimated  using  the  results  for  random 
walk  in  two  dimensions,  <r2>  =  4Dt,  where  r  is  a  characteristic  length,  in  this 
case  the  beam  radius(~3  pm),  and  t  is  the  decay  time  of  the  observed  intensity 
fluctuations. 

The  SH  intensity  was  recorded  every  0.5  s  for  half  an  hour,  p-nitrophenol 
solution  is  used  as  the  subphase  as  it  is  a  good  probe  of  the  surface  generating 
large  SH  signal.  Experiments  were  performed  at  37,  50,  100,  200  and  400 
A2/molecule  of  Cl  6  aniline.  The  results  are  shown  in  figure  3. 5, next  page. 
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The  fluctuations  in  the  SH  signal  from  400A2  and  37A2/molecule  is  different 
from  that  generated 

by  the  Cl  6  aniline  at  coverages  of  400,  200  and  100A2/molecule.  This  can  be 
understood  in  terms  of  the  phase  diagram.  At  37A2/molecule  the  Cl 6  aniline  on 
the  surface  is  in  a  homogeneous  liquid-like  phase,  while  the  other  three 
measurements  correspond  to  its  being  inside  the  gas-liquid  coexistence  region. 
At  very  low  coverages  ('gas'),  the  SH  signal  is  the  same  as  for  a  bare  water 
surface.  However,  as  the  Cl 6  aniline  packing  is  increased,  a  substantially 
larger  SH  signal  is  generated.  It  is  notable  that  the  'low'  signal  in  200A2  :s  tho 
same  signal  as  400A2.  This  signal  corresponds  to  patches  of  the  surface  with 
the  low  density  gas-like  coverage  in  the  irradiated  region.  The  'high'  signals,  on 
the  other  hand,  correspond  to  the  high  density  'liquid'-like  coverage  in  the 
irradiated  region.  At  200A2/molecule  the  coverage  is  low,  and  there  are  few 
liquid-like  islands  dispersed  in  gas-like  coverage.  As  the  packing  is  increased, 
at  100  and  50A2/molecule,  more  'high'  signals  are  obtained,  corresponding  to  a 
higher  fraction  of  the  Cl  6  aniline  in  a  liquid-like  state.  Finally,  at  37A2/molecule, 
the  surface  is  now  covered  by  a  homogeneous  layer  of  palmitic  acid  in  its 
'liquid-like'  phase.  The  signal  fluctuations  for  both  the  gas  phase  surface  400A2 
and  the  homogeneous  liquid  37A2  are  ascribed  to  random  noise  as  discussed 
below. 

As  a  further  test  of  the  significance  of  the  observed  fluctuations,  the 
probability  density  function  for  our  data  was  plotted  in  figure  3.6.  The  probability 
density  functions  of  400  and  37A2  have  been  fit  to  a  poisson  distribution  which 
shows  that  fluctuations  observed  in  these  cases  is  random  noise,  arising  from  a 
homogeneous  surface.  However  the  distributions  in  200,  100,  50A2do  not  look 
like  simple  ones  and  can  probably  be  fitted  to  a  linear  combination  of  poissons, 
which  is  consistent  with  the  notion  of  a  heterogeneous  surface. 

The  autocorrelation  <A1(0)  Al(t)  >  vs.  time  are  shown  in  figure  3.7.  The 
data  for  50,  100  and  200  A2/molecule  show,  a  monotonic  decay  in  the 
correlation  function  while  those  for  400  and  37A2  show  only  a  delta  function 
indicating  a  randomly  uncorrelated  noise. 
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Since  the  high  signals  obtained  for  50,  100,  and  200A2/molecule,  i.e. ,  in 
the  coexistence  reqion,  are  of  the  same  magnitude  as  the  signal  obtained  for 
the  pure  liquid  state  e.g.  at  37A2/molecu!e,  we  conclude  that  the  clusters  are 
comparable  in  size  with  our  incident  beam.  By  fitting  our  data  to  a  diffusion 
model,  a  diffusion  constant  of  ~  10'8  cm2/s  and  a  size  of  ~1pm  are  obtained. 

3.4  Orientation  Phase  Transition  Studied  by  SH  Fluctuation 

We  now  wish  to  report  the  observation  of  orientational  fluctuations  of  lipid 
molecules  at  the  air/water  interface,  which  we  attribute  to  a  heretofore  unknown 
orientational  phase  transition.  These  fluctuation  were  shown  to  be  orientational 
by  the  observations  that  the  they  depended  on  the  polarization  of  the  incident 
light  field  and  the  polarization  of  the  second  harmonic  light  field  that  was 
detected.  The  experimental  results  of  this  study  bear  on  the  organisation  of  lipid 
molecules  at  the  air/water  interface,  a  subject  that  has  been  and  remains  one  of 
intense  experimental  and  theoretical  activity.  The  molecule  we  studied  is 
hexadecylaniline,  a  long  chain  aniline  where  the  aniline  head  group  has  C2v 
symmetry.  For  an  interface  that  is  isotropic  in  the  surface  plane,  i.e.  is 
rotationally  invariant  with  respect  to  the  interface  normal,  populated  by 
molecules  of  C2v  symmetry  the  explicit  relation  between  x(2)  and  a(2)  yields 
for  two  of  the  x^2)  elements  that  we  will  be  concerned  with  the  following 
expressions, 

x zxx=  2  Ns[f<aZzr^xx’axzx'e''l')  +  °2xx<cos8>)  <eq  3  7a) 

Xxzx=  2  Ns[f(azzz’'4xx-C4z>x’e’'t')  +  <\zx<cos9>1  <eq  3  7b> 

where  Ns  is  the  lipid  surface  density  (molecules/cm2)  and 

f^azzz’azxx’axzx’e-'l')  =  <*jzz<cos9sin2e>  •  (“j5,x+2'\zx><cos95in20sin2V>. 
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6  is  the  angle  between  the  molecular  symmetry  axis  z  and  the  surface  normal  Z, 

and  v  is  the  rotation  angle  of  the  molecule  about  its  symmetry  axis  z.  From  Eqs. 

(2)  (2) 

3.7  and  3.8  we  see  that  \zx  and  Xzxx  differ  only  in  the  nonlinear  polarizability 
coefficients  of  <cos8>. 

(2)  (2) 

The  SH  signal  corresponding  to  these  two  elements,  \zx  and  Xzxx 

(Eqs  3.7  and  3.8),  were  measured  as  a  function  of  time.  On  further  compression 

to  densities  above  37  A2  where  the  monolayer  is  in  the  homogeneous  liquid 

(2) 

phase  region,  a  non-Poisson  fluctuations  were  observed  when  Xxzx  is 

(2) 

measured.  However,  for  Xzxx  element,  the  fluctuations  observed  were  very 

close  to  a  Poisson  distribution.  We  therefore  conclude  that  we  are  detecting  a 
polarization  dependent  and  hence  an  orientational  fluctuation.  It  should  be 
noted  that  there  is  no  indication  from  the  k-A  phase  diagram  of  a  phase 
transition  in  this  density  region  of  the  liquid  phase.  On  compressing  the 
monolayer  further  we  find  that  fluctuations  decrease  and  at  28A2  the 
fluctuations  are  approximately  Poisson  in  character.  We  were  unable  to  go  to 
higher  coverages  because  of  the  collapse  of  the  monolayer  below  28A2. 

We  now  consider  the  basis  for  the  orientational  fluctuations.  Since  the 
SH  signal  is  due  chiefly  to  aniline  we  note  that  the  fluctuations  that  we  observe 
must  be  due  to  the  fluctuations  in  the  orientation  of  the  aniline  head  group. 
However  one  cannot  conclude  from  this  that  the  driving  force  for  these 
fluctuations  are  due  only  to  head  group  interactions,  e.g.  via  a  dipole-dipole 
coupling  between  the  head  groups.  It  is  quite  reasonable  that  chain-chain 
interactions  at  these  relatively  high  densities  play  some  role  in  driving  the  head 
group  orientational  phase  transition.  We  know  that  the  chain-chain  interactions 
cannot  be  the  sole  factor  however  since  we  do  not  observe  any  orientational 
fluctuations  for  an  equal  chain  length  lipid  with  a  different  head  group,  namely 
hexadecylphenol  (n-Ci6H33-C6HsOH). 

Let  us  now  consider  the  orientation  of  the  aniline  head  group,  which  is 
defined  by  0,  the  angle  of  the  aniline  symmetry  axis  with  the  interface  normal  Z 
and  v/  the  angle  of  twist  about  the  aniline  symmetry  axis.  Noting  that  \j r  is 
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present  only  in  the  function  f°r  *xzx  *zxx’  ^qs. 


3.7  and  3.8,  it  appears  that  the  0  dependence  in  the  aXZx<cos0>  term  is 


(2) 


responsible  for  the  fluctuations,  since  it  is  only  this  term  that  differs  for  X^zx  and 


,<2> 

LZXX- 


.(2) 


(2) 


From  theoretical  estimates  axzx  is  much  larger  than  a7YY  for  aniline. 


zxx 


This  would  explain  why  we  readily  observed  orientational  fluctuations  in  the 

(2)  (2) 

Xxzx  element  and  not  in  the  >^xx  element. 


In  order  to  examine  the  spatial  range  of  the  orientational  fluctuation  we 
varied  the  laser  beam  area  by  a  factor  of  four  from  a  radius  of  3pm  to  6pm.  The 
results  of  autocorrelation  times  for  the  two  beam  areas  for  various  monolayer 
coverages  is  given  in  Table  1. 


36A2 

34A2 

32A2 

30A2 

28  A2 

3pm 

9 

1  3 

16 

1  7 

~1 6 

6pm 

5 

8 

1  2 

1  7 

~15 

Table  I. 

Correlation  Time  tc(sec)  Dependence  of  Area  per  Molecule 
and  Observation  Spot  Size 


As  the  laser  spot  size  is  increased  the  number  of  domains  viewed  increases, 
and  therefore  a  more  rapid  decorrelation  would  be  expected.  As  the  transition 
point  is  approached  the  domains  increase  in  size  and  the  decay  of  the 
autocorrelation  function  would  be  slower,  though  still  showing  a  change  with 
spot  size.  At  30A2  the  system  is  near  the  transition  point  and  the  domains  are 
very  large,  becoming  larger  than  the  beam  size.  This  is  seen  in  the  result  that 
the  same  decay  time  is  observed  for  the  3p  and  6p  radii  beam,  which  indicates 
that  the  orientational  correlation  length  is  greater  than  6p.  If  the  transition  was 
second  order  then  this  correlation  length  would  diverge  at  the  transition  point. 
Based  on  these  scaling  experiments,  which  indicate  a  long  correlation  length, 
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we  infer  that  the  phase  transition  is  likely  to  be  weak  first  order  if  not  second 
order. 

As  alluded  to  earlier  the  importance  of  chain-chain  interactions  on  the 
observed  orientational  fluctuations  is  seen  in  the  strong  dependence  of  the  SH 

fluctuations  on  the  chain  length  of  the  aniline  lipid.  The  autocorrelation  functions 

(2) 

and  probability  density  functions  of  \zx  for  the  Cl 2,  Cl 4,  Cl  6  anilines  at  a 

coverage  of  34A2  are  shown  in  figure  3.8.  It  is  evident  from  figure  5  that  the 
shorter  chain  anilines  have  shorter  correlation  times  and  distributions  closer  to 
Poisson,  thus  indicating  a  more  homogeneous  surface  on  the  time  scale  of  our 
measurements.  In  fact  Cl  2-aniline  does  not  show  any  deviation  from  Poisson 
fluctuation  for  the  experimental  sampling  time  scale  of  0.5  sec  that  we  use.  This 
could  be  due  to  a  correlation  length  that  is  much  smaller  than  the  laser  beam 
area  (3|i  radii),  and/or  a  relaxation  that  is  much  faster  than  the  0.5  sec  sampling 
time.  It  is  perhaps  not  surprising  that  the  stronger  attractive  interactions  for  the 
longer  chain  anilines  would  be  more  restrictive  of  the  head  group  orientational 
motions. 

To  examine  the  effect  of  the  head  group  on  these  orientational 
fluctuations  we  studied  Cl  4,  Cl  6  and  Cl  8-phenols  as  well  as  Cl  6  anilinium 
and  Cl 8  phenolate  syetem.  Although  phenol  and  aniline  have  similar  dipole 
moments,  and  the  chain  lengths  were  the  same  for  both  the  aniline  and  phenol 
lipids  we  observed  no  orientational  fluctuations  for  any  of  the  phenols.  The 
observation  indicates  that  the  interactions  of  the  head  groups  with  the 
surrounding  water  molecules  is  an  important  factor,  recognizing  that  the 
hydrogen  bonding  properties  of  the  aniline  vs  phenol  are  different.  The 
necessary  balance  between  the  various  forces,  chain-chain,  hydrogen  bonding, 
and  electrostatic,  involved  in  the  phase  transition  is  demonstrated  by  our 
observation  that  the  charged  lipid  C18-phenolate  (CI8-C6H5O')  manifests 
orientational  fluctuations  similar  to  the  Cl  6-aniline.  The  charged  anilinium  lipids 
(Cl  6-C6H5NH3+)  and  as  noted  earlier  the  neutral  phenol  lipids  do  not  show 
orientational  fluctuations.  It  is  of  interest  that  the  neutral  and  charged  lipids 
showing  these  fluctuations  have  similar  k- A  phase  diagrams. 
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4.  Solid-Liquid  and  Electrified  Interfaces 

4.1  Studies  of  Silica/electrolyte  solution  interface  by  SHG  The 

glass-water  interface  is  an  important  example  of  a  solid-liquid  interface.  In  the 
case  of  silica  in  contact  with  electrolyte  solution,  such  as  NaCI,  the  surface 
complexation  dominates  at  higher  electrolyte  concentration  [12,13]  .  Thus,  the 
main  surface  reaction  at  silica/electrolyte  (NaCI)  solution  is 

Kint 

SiOH  +  Nas+  SiQ-Na++  Hs+  (eq4.1) 


With  the  intrinsic  surface  equilibrium  constant 

K'NT  “  tSiONa  3hT  eXPt'e('*V  (eq  4  21 


Where  SiONa+  is  an  ion  pair  complex,  simplified  as  SiONa,  while  fsiOH, 
fsio-and  fsiONa  are  the  fraction  of  SiOH,  SiO'  and  SiONa,  respectively,  an+  and 
aNa+  are  activities  of  proton  and  sodium  ion,  respectively.  4^  is  the  electrostatic 
potential  affecting  the  SiO'  charge  sites, which  are  assumed  to  be  located  in  the 
same  zero  plane  as  that  for  H+  and  OH'.  Invoking  the  Stern  model  of  an 
electrified  interface,  the  role  of  counterions  is  included  in  the  potential  4/p  . 

The  basic  features  of  the  SH  signal  from  silica/electrolyte  interface  are 
that  at  fixed  ionic  strength,  the  SH  signal  increases  with  increasing  pH  ,  while  at 
fixed  pH,  SH  signal  decreases  with  increasing  concentration  of  NaCI  .  It  is  well 
known  that  surface  silanols  dissociate  more  and  more  with  increasing  pH  and 
that  there  is  more  surface  complexation  with  increasing  concentration  of  NaCI. 
Thus,  we  may  conclude  that  the  SH  signal  from  silica/electrolyte  interface 
increases  with  increasing  total  surface  charge  density  (Nsio-  +NsiONa  )  and  that 
at  the  same  surface  charge  density  ,  the  SH  signal  decreases  with  more  surface 
complexation.  The  SH  signal  may  arise  from  two  different  mechanisms,  either 
from  the  second  order  hyperpolarizability  or  the  third  order  hyperpolarizability 
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arising  from  the  electrical  potential  across  the  solid/liquid  interface  due  to  the 
surface  charge  [1,14],  or  combination  of  both. 

At  higher  electrolyte  concentration,  surface  complexation  dominates  and 
the  amount  of  SiONa  is  much  larger  than  that  of  SiO-,  e.  g.  f  sioNa  » fsiO- 
Under  this  condition,  the  two  different  mechanisms  end  up  with  the  identical 
relationship  between  the  SH  signal  and  the  density  of  the  surface  charge 
species,  e.g.  SiO'  and  SiONa.  It  was  found  that  at  a  given  pH,  the  fraction  of 
SiONa,  fsiONa  is 


fSiONa  = 


(lS)1/2-  ds  (PH  =  2))1/2 
0s(pH  =  13))1/2  -  (Is  (pH  =  2))1/2 


(eq  4.3) 


Where  (ls  )1/2,  (ls  (Ph=2))1/2  and  (lS(pH=i3))1/2  are  the  square  roots  of  the 
SHG  intensity,  respectively.  This  is  not  surprising  since  the  surface 
complexation  dominates  and  the  net  charge  density  in  diffuse  layer  is  very  small 
compared  to  the  charge  density  at  the  zero  plane,  which  is  linearly  proportional 
to  E0,  the  electric  field  at  the  interface.  It  is  thus  fortunately  that  we  can  obtain 
unique  results  of  equilibrium  constants  of  surface  reaction  (eq.  4.1)  by  fitting  the 
experimental  data  even  if  we  could  not  separate  these  two  different 
mechanisms,  although  it  is  of  importance  and  interest  to  separate  them.  An 
attempt  to  separate  these  two  mechanisms  by  varying  the  temperature  is  in 
progress. 

Knowing  the  fraction  of  the  ion-pair  complexes  SiONa,  fsiONa  from  the 
experimental  measurement  according  to  eq.  4.3  ,  we  can  determine  the  intrinsic 
surface  complexation  constant  Kint  of  reaction  4.1 .  The  plot  of  fsiONa  vs  bulk 
pH  looks  like  a  titration  curve  with  two  pKa  values.  This  suggests  that  there  are 
two  different  exchange  sites  with  distinct  intrinsic  equilibrium  constants  K|NT. 
Thus,  a  two-site  ion  exchange  model  proposed  by  Allen  et  al  [15]  was  adopted 
to  analyze  our  data.  The  two-site  model  assumes  that  two  distinct  types  of 
silanol  groups  (I  and  II)  are  available,  and  that  the  intrinsic  surface  complexation 
constants  K|Nj  differ  so  much  that  the  type  sites  I  first  filled  completely,  and  that 
only  then  do  the  type  II  sites  fill. 

The  result  is  shown  in  following  figure  4.1 .  It  was  found  that  about  20% 
of  the  surface  silanols  have  a  pKinj  =  5.3,  while  the  other  80%  of  silanols  with 
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PKint  =  9-2.  This  is  in  excellent  agreement  with  the  results  obtained  by  Allen  et 
al,  with  silica  gel  solution  by  potential  titration  method  [15]  . 


Figure  4.1 

Fraction  of  SIONa  versus  bulk  pH 
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4.2  Studies  of  Molecular  Properties  at  the  Surface  of  a  Liquid  Jet 
by  Second  Harmonic  Generation  We  report  here  on  the  use  of  surface 
second  harmonic  generation  to  study  molecular  equilibrium  and  dynamic 
processes  at  the  surface  of  a  liquid  water  jet. [16]  The  use  of  a  liquid  jet  is  a 
potentially  important  way  to  study  interfaces  in  that  it  permits  the  study  of 
equilibrium,  e.g.,  orientational  structure  and  solute  adsorption  energetics,  as 
well  as  dynamic  molecular  processes  in  systems  where  the  interfering  effects  of 
irreversible  photochemical  and  thermal  processes  can  be  minimized  or  avoided 
by  the  continuously  generated  fresh  surface.  It  also  offers  the  possibility  of 
generating  dpsired  nonequilibrium  populations  at  an  air/liquid  interface  by 
controlling  the  characteristics  of  the  solid/liquid  interface  of  the  nozzle 
immediately  prior  to  the  formation  of  the  jet.  In  this  study,  it  is  the  dynamics  of  the 
adsorption  of  an  organic  molecule,  p-nitrophenol,  to  the  air/water  jet  interface 
that  we  report. 

The  initial  population  density  of  nitrophenol  molecules  at  the  liquid  inter¬ 
face  that  exit  the  glass  nozzle  is  determined  by  the  nitrophenol  interactions  at 
the  glass/water  interface,  not  its  interactions  at  the  later  developing  air/water  jet 
interface.  Thus,  the  air/water  jet  interface  has  a  surface  density  of  nitrophenol 
molecules  that  initially  is  not  its  equilibrium  value.  The  experiments  measure  the 
time  it  takes  the  interface  to  achieve  its  equilibrium  physical  and  chemical  state. 
From  our  determination  of  the  jet  velocity  we  know  the  elapsed  time  from  the  jet 
nozzle  to  each  position  along  the  axis  of  the  jet.  The  time  dependent  change  in 
the  nitrophenol  population  at  the  air/water  jet  interface  can  therefore  be  ob¬ 
tained  by  measurement  of  the  SH  signals  at  various  distances  from  the  jet  noz¬ 
zle.  We  take  the  square  root  of  the  SH  signal,  since  for  a  coherent  process  this 
is  the  quantity  that  is  proportional  to  the  number  of  nitrophenol  molecules/cm2  at 
the  interface.  Using  this  technique  we  have  been  able  to  study  the  kinetics  and 
energetics  of  adsorption  of  p-nitrophenol  to  the  air/water  interface.  We  found 
that  the  adsorption  kinetics  to  the  liquid  jet/air  interface  could  be  described  by  a 
dynamic  Langmuir  model.  From  the  kinetics  we  obtained  the  rate  constants  of 
4.4  ±  0.2  x  104  sec1  for  adsorption  and  6  ±  2  sec1  for  desorption.  The  free  en¬ 
ergy  of  adsorption  for  the  process  was  found  to  be  AG°ads  =  -  5.1  ±  0.2 
Kcal/mole,  which  is  the  same  as  for  the  static  (equilibrium)  case  .  The  kinetic 
data  also  yields  the  initial  surface  population  of  p-nitrophenol  at  the  air/water 
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interface,  which  thereby  gives  information  about  the  adsorption  at  the  preceding 
glass/water  interface  of  the  jet  nozzle.  We  also  found  that  adsorption  of  p- 
nitrophenol  to  the  air/water  interface  is  considerably  slower  than  diffusion 
controlled,  thus  indicating  that  the  adsorption  rate  is  kinetically  controlled  by  a 
barrier.  The  orientation  of  nitrophenol  at  the  liquid  jet  /air  interface  was  the  same 
as  for  the  static  (equilibrium)  liquid/air  interface,  which  indicates  that 
orientational  equilibrium  was  achieved  prior  to  the  earliest  measurement  time, 
0.55  msec.  In  addition,  we  note  that  this  method  offers  the  possibility  of  generat¬ 
ing  various  nonequilibrium  populations  of  species  at  air/liquid  interfaces  and 
also  opens  the  possibility  of  investigating  static  and  dynamic  interfacial  pro¬ 
cesses  by  avoiding  the  interfering  effects  of  irreversible  photochemical  and 
thermal  processes  common  to  many  systems. 

4.3  Phenols  and  Anilines  at  the  Air/Water  interface  Surface  second 
harmonic  generation  has  been  used  to  study  the  adsorption  of  alkyl  phenols 
and  anilines  and  their  respective  ions  to  the  air/water  interface. [17]  The  free 
energies  of  adsorption  of  the  neutral  molecules  were  obtained  using  a  simple 
Langmuir  model.  The  quantitative  effects  of  the  opposing  hydrophobic  and 
solvation  forces  on  molecular  adsorption  were  investigated  by  varying  the  chain 
length  of  alkylphenolate  and  alkylanilinium  ions.  An  alkyl  chain  of  at  least  five 
carbons  was  found  to  be  necessary  to  balance  the  bulk  ionic  solvation  and 
thereby  bring  the  ion  interface  population  to  a  level  above  the  water 
background.  The  adsorption  of  the  ions  can  be  quantitatively  described  by 
using  a  modified  Langmuir  model  that  includes  the  electrostatic  free  energy  of 
the  electric  double  layer  as  given  by  the  Gouy-Chapmann  theory.  From  the 
excellent  fit  of  the  second  harmonic  data  to  this  treatment,  the  chemical  and 
electrostatic  free  energies  of  adsorption,  interface  electric  potential,  and  in¬ 
terface  pH  were  obtained. 
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5.  Ultra-fast  Dynamics  at  the  Liquid/Air  Interface  Studied  by 
Second  Harmonic  Generation  Techniques. 

5.1  Photoisomerization  Dynamics  as  a  Probe  for  the  Interface  The 

different  dielectric  (polarity)  and  transport  properties  (viscosity)  can  result  in  very 
different  static  and  dynamic  properties  in  the  bulk  and  interfacial  regions.  To 
probe  the  liquid/vapor  interface,  we  used  the  surface  selective  technique  of 
second  harmonic  generation  (SHG).  We  determined  how  chemical  reactions 
differed  in  the  air/water  interface  compared  with  the  bulk.  We  carried  out 
picosecond  time  resolved  measurements  of  excited  state  lifetimes  of  reacting 
molecules  in  the  air/water  interface.  In  these  pump/probe  studies  we  used 
second  harmonic  generation  methods  to  probe  the  dynamics  in  the  interface. 
We  also  measured  the  bulk  dynamics  of  the  same  process  in  the  same  solution 
using  picosecond  time  resolved  fluorescence.  In  this  way  we  obtained 
information  not  only  about  the  interface,  by  determining  how  it  differs  from  bulk 
dynamics,  but  also  about  the  molecular  process  itself. 

We  earlier  reported  on  the  photoisomerization  of  the  cyanine  dye  DODCI 
in  the  water/air  interface. [18]  We  measured  the  SH  signal  as  a  function  of  time 
delay  between  the  pump  and  probe  pulses.  From  this  study  we  determined  that 
the  photo-induced  isomerization  time  constant  is  220  ±  40  ps  in  the  air/water 
interface  versus  520  ±  60  ps  for  the  bulk  region  of  the  same  solution.  The  faster 
rate  at  the  interface  could  indicate  a  smaller  friction  for  motion  along  the 
isomerization  coordinates  in  the  interface. 

In  addition  to  these  important  kinetic  results,  we  also  found  that  the  ratio 
of  the  yield  of  trans-photoisomer  compared  to  the  cis-isomer  was  dt/c  =  40-50%, 
which  was  the  same  as  that  found  in  the  bulk  aqueous  phase.  Obviously,  the 
ratio  of  product  yield,  i.e.,  [cis-isomer]/[trans-isomer],  does  not  depend  on  the 
barrier  crossing  dynamics  of  the  excited  Si  to  twisted  intermediate  state,  but 
rather  depends  directly  on  the  way  the  twisted  polar  intermediate  decays.  The 
decay  of  the  twisted  polar  intermediate  into  either  the  cis-form  or  the  trans-form 
depends  on  the  solvent  environment.  Based  on  our  dynamic  measurements  it 
can  be  argued  that  the  two  ground  state  isomers  in  the  interfacial  region 
experience  a  similar  solvation  environment,  not  unlike  that  as  found  in  the  bulk 
aqueous  phase. 
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Figure  5.1 

Photoisomerization  of  DODCI  (3,3'-dlethyloxadlcarbocyanlne  Iodide) 


To  further  explore  the  influence  of  the  chemical  composition  of  the  air- 
liquid  interface  on  photoisomerization  reactions  we  now  have  extended  our 
studies  of  the  cyanine  dye  DODCI  at  the  air/water  interface  to  include  the  effect 
of  nonpolar  molecules  at  the  interface.  By  varying  the  chemical  composition  of 
the  air/liquid  interface  a  corresponding  change  in  the  photophysics  of  the 
adsorbed  cyanine  dye  molecule  was  seen.  Specifically,  we  measured  the 
photoisomerization  of  DODCI  at  the  surface  of  a  spread  monolayer  of  hexadecyl 
alcohol.  The  lipid  surface  provides  DODCI  with  a  nonpolar  environment.  It  was 
found  that  the  surface  adsorbed  DODCI  molecule  undergoes  rapid 
isomerization  (=  300  ps)  and  <|>t/c  <  10%.  The  fast  kinetics  for  the  air-lipid 
interface  is  similar  to  that  found  for  the  air-water  interface,  but  is  faster  than  than 
the  =  1  nsec  time  expected  for  a  nonpolar  bulk  liquid  phase.  The  similar  air- 
liquid  photoisomerization  kinetics  supports,  however,  our  earlier  suggestion  that 
the  interfacial  friction  is  inherently  low  compared  to  what  is  found  in  the  bulk 
liquid  phase.  Finally,  the  dramatic  reduction  in  <f>t/c  for  the  lipid  rich  interface 
(relative  to  the  water  surface)  is  consistent  with  the  solvation  properties  of  the 
lipid  for  the  DODCI  molecule.  It  is  known  that  0t/c  in  long  chain  alcoholic 
solvents  is  exceedingly  low  (<  10%),  just  as  what  we  determined  for  the  air-lipid 
interface. 
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5.2  Dynamics  of  Proton  Transfer  at  the  Interface  We  also  investigated 
the  dynamics  of  proton  transfer  reactions  at  the  air/water  and  lipid/water 
interface.  At  present  there  is  little  direct  information  on  this  fundamental  reaction 
in  the  interfacial  region.  We  made  direct  kinetic  measurements  using  the  pump- 
probe  SHG  method  already  outlined.  Specifically,  we  examined  the  fast  proton 
transfer  reaction  of  aromatic  phenols  (such  as  b-naphthol  and  p-nitrophenol)  as 
well  as  aromatic  amino  compounds  (aminopyrene).  From  bulk  work  we  know 
that  in  the  excited  state  these  molecules  are  more  acidic  (pKa*  <  <  pKa).  By 
photoexciting  the  compounds  at  the  interface,  we  induced  an  acid-base 
reaction,  which  was  detected  by  the  time-dependent  change  in  the  SH  signal. 

From  the  recovery  kinetics  following  photobleaching  of  aminopyrene  at 
low  pH  (<  2)  we  conclude  that  interfacial  deprotonation  to  form  the  neutral 
molecule  is  rapic  (<  320  ps) 


Figure  5.2 

proton  transfer  in  aminopyrene 


ps).  The  deprotonation  kinetics  in  the  air/water  interfacial  region  are  similar  to 
what  is  known  in  the  bulk  aqueous  phase.  Another  example  of  a  proton  transfer 
reaction  involves  the  deprotonation  of  a  neutral  molecule  to  form  the  charged 
anion,  as  exemplified  by  p-nitrophenol.  Again  we  found  that  deprotonation  time 
is  extremely  rapid,  taking  place  in  less  than  50  ps.  This  result  is  important,  since 
we  earlier  established  that  under  equilibrium 
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Figure  5.3 
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Proton  transfer  In  PNP 

conditions  the  anion  is  repelled  from  the  surface,  i.e.,  it  does  not  adsorb  to  the 
surface  but  remains  in  the  bulk  water.  The  fact  that  excited  state  deprotonation 
(to  produce  the  charged  species)  does  in  fact  occur  at  the  surface  sheds 
important  light  on  the  interfacial  solvation  of  the  solute  molecule.  Presumably 
the  phenolic  chromophore  is  directed  sufficiently  toward  the  bulk  water 
molecules  so  that  they  can  readily  accept  the  proton  from  the  photoexcited 
phenol,  and  thereby  achieve  charge  separation. 

5.3  Rotational  Motion  at  the  Air/Water  Interface  From  SHG 

experiments,  we  measured  specific  elements  of  the  surface  nonlinear 
susceptibility,  %ijk  Figure  5.4  shows  the  time  dependence  of  xzxx  and  xxzx-  At 

negative  time  (t  <  0),  before  the  pump  pulse  has  arrived,  we  are  probing  the 
ground  state  molecules.  Due  to  photoexcitation,  molecules  change  their 
nonlinear  polarizability,  thus  the  SH  signal  changes.  In  the  present  case,  we 
detected  a  signal  decrease  for  both  elements,  which  we  usually  call  bleaching. 
After  the  initial  decrease  in  signal  due  to  photoexcitation,  xzxx  increases  and  at 

sometime  t  it  goes  even  beyond  the  value  corresponding  to  negative  time. 
Then  the  signal  decays  back  to  the  initial  value.  However  for  xxzx.  after  the 

initial  decrease  due  to  photoexcitation,  there  is  a  further  decrease  in  its  signal, 
which  later  also  decays  back  to  the  initial  value.  As  it  can  be  seen,  our  data 
clearly  show  a  polarization  dependence  of  the  dynamics.  Different  elements  of 
the  SH  signal  show  a  different  time  dependence  for  the  relaxation  process. 
This  indicates  that  rotational  motions  have  taken  place. [19]  If  there  were  no 
rotational  motions  involved,  the  time  dependent  SH  signal  for  xzxx  and  xxzx 

would  only  come  from  the  population  term,  so  the  dynamics  curve  should  be  the 
same  for  both  elements,  just  a  population  relaxation.  However  this  is  clearly  not 
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the  case.  Since  the  SH  signal  is  polarization  dependent,  we  believe  that  we 
are  observing  rotational  motions  at  the  air/water  interface. 

Now,  how  can  we  explain  the  initial  decrease  in  signal  due  to 
photoexcitation  that  we  have  seen  in  both  elements?  At  negative  time,  before 
the  pump  pulse  has  arrived,  we  are  probing  ground  state  molecules,  Xjjk(t  <  0). 

Right  after  the  excitation,  but  before  any  population  or  orientational  relaxation 
has  taken  place,  SH  signal  can  be  expressed  as:  xjjk(t=0), 

Xijk(t  <  0)  =  N0  <  a(2)  > 

Xjjk(t  =  0)  =  No  <  a(2)  >  -  Ns*  (<  a(2)  >  -  <  a  >) 

(eq  5.1) 

From  our  experimental  data,  the  SH  signal  decreased  due  to  photoexcitation. 
This  implies  that  the  nonlinear  polarizability  of  the  excited  state  molecules  is 
smaller  than  that  of  the  ground  state  molecules. 

There  is  another  interesting  feature  in  our  data:  at  sometime  t,  the  SH 
signal  for  xzxx  goes  even  beyond  the  initial  value  at  negative  time.  The 
difference  in  signal  for  xzxx  at  time  t  and  at  negative  time  can  be  expressed  as: 

AXZxx  =  Ns*(t)  (<  a  (t)>*  a  (2)  -  <  Q  (t)>  cc(2> )  (eq  5.2) 

where  <  Q  (t)>  is  the  orientational  factor  which  is  related  to  the  direction  cosines 
in  terms  of  the  orientational  angles.  Our  data  show  that  AXZxx>  0  at  time  and 
we  know  that  the  nonlinear  polarizability  of  the  excited  state  is  smaller  than  that 
of  the  ground  state.  Thus  the  orientational  factor  for  the  excited  state  molecules 
has  to  be  larger  than  that  for  the  ground  state  molecules.  This  implies  that  the 
ground  and  excited  states  molecules  have  different  final  orientational 
distributions.  A  similar  argument  can  be  applied  to  explain  the  further  decrease 
in  xxzx  after  the  bleaching.  Actually,  if  we  substitute  the  direction  cosines  for  a 
set  of  Euler  angles,  the  orientational  factor  indeed  increases  for  xzxx  and 
decreases  for  xxzx- 

In  summary,  the  features  of  our  data  are:  there  are  two  parts  in  the  kinetic 
curves.  A  polarization  dependent  part  at  early  times,  where  the  SH  signal 
increased  in  xZxx  dut  decreased  in  xxzx  We  think  this  part  is  due  to  the 

rotational  motions  at  the  surface.  As  it  can  be  seen,  this  takes  place  between 
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500  to  1000  picoseconds.  A  slow  decay  which  is  similar  in  both  elements  is 
probably  due  to  the  excitation  relaxation  of  Rh6G  molecules  at  the  surface.  It  is 
in  the  time  scale  of  ~  4  ns  which  is  comparable  to  the  lifetime  of  Rh6G  in  the 
bulk  [20]. 

6.  Femtosecond  Laser  Studies  of  Electrons  in  Liquids 

6.1  Introduction  The  properties  of  electrons  in  liquids,  ranging  from 
hydorcarbons  to  water,  are  of  widespread  interest  and  are  recogized  to  be 
important  in  a  wide  variety  of  physical  problems  such  as  electrons  in 
amorphous  solids,  liquid  state  chemistry,  and  biological  problems  such  as 
electron  transfer.  We  have  been  uisng  femtosecond  UV  laser  pulses  to  directly 
study  the  dynamics  of  geminate  recombination  of  electrons  and  cations  and 
the  solvation  dynamics  in  neat  water.  The  electrons  were  generated  by 
mulitphoton  ionization  of  neat  alkanes  and  neat  water  was  as  well  as  by 
multiphoton  photodetachment  of  electrons  from  halide  ions  and  the  hydroxide 
ion  in  aqueous  solutions. 

6.2  Electron  Solvation  Dynamics  in  Neat  Water  There  are  two  contrasting 
descriptions  of  electron  solvation  dynamics  that  have  been  proposed,  the  dielectric 
continuum  model  (DCM)  and  the  two-state  model.  The  dielectric  continuum  model 
(DCM)  ,  has  been  used  extensively  in  simple  analytic  theories  of  various  chemical 
processes  in  polar  liquids.  [21]  It  predicts  that  a  continous  absorption  spectra  blue 
shift  corresponding  to  the  solvent  molecular  rotation.  The  two  state  model  means  that  a 
wet  electron,  identifiable  by  its  characteristic  absorption  spectra,  undergoes  a 
transition  to  the  solvated  electron.  In  other  words,  there  are  only  two  electron  species 
present  in  the  kinetics  rather  than  a  series  of  wet  electron  species,  having  their  own 
spectra,  that  evolve  into  the  equilibrium  solvated  electron.  In  the  two  state  model,  if  the 
absorption  of  the  two  species  (states)  overlap  then  there  must  be  a  wavelength  where 
the  two  spectra  intersect  and  therefore  have  the  same  absorption  coefficient  e.  A 
consequence  of  the  presence  of  an  isosbestic  point  can  be  seen  by  considering  the 
absorption  A(t)  at  a  time  t  and  a  wavelength  X. 


(eq.  6.1) 

3  0 


A(X,t)=£wel(X)Nwe,(t)  +  EsoA)Nsoi(t) 
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where  Nwet(t)  and  NS0|(t)  are  the  concentrations  of  wet  and  solvated 
electrons.  At  the  isosbestic  wavelength,  £Wet=£soi  and  the  kinetics  will  therefore  be 
constant  after  the  wet  electron  is  generated  since  the  total  electron  concentration, 
Nwet(t)  +  NSOi(t),  is  conserved. 

To  find  the  isosbestic  point  we  scanned  from  the  near  infrared  to  the  visible.  The 
time  evolution  of  the  aqueous  electron  at  several  different  probe  wavelengths  is  shown 
in  Figure  6.1,  see  next  page.  For  absorption  signal  at  625  nm,  the  monotonically  rising 
component  corresponds  to  the  appearance  of  the  wet  or  pre-solvated  electron  and  its 
transition  into  the  fully  solvated  electron.  At  1000  nm,  the  absorption  increases  quickly 
then  decays  more  slowly.  We  attribute  this  decay  to  the  transition  from  the  wet  electron 
to  the  solvated  electron.  At  820  nm,  there  is  a  rise  in  the  absorption,  corresponding  to 
the  wet  electron  formation,  then  a  distinct  leveling  off  of  the  kinetics.  We  observed 
the  isosbestic  wavelength  kinetic  behaviour  at  820  nm,  thus  we  concluded  that  the 
wet  to  solvated  electron  transition  involves  only  two  states. 

Accordingly,  we  used  a  simple  kinetic  model  to  fit  the  data. 


ki  k2 
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(eq.  6.2) 


The  quasi-free  electron  absorption  is  neglected  due  to  its  extended  nlane 
waves  characterization.  Using  the  known  values  for  the  solvated  electron  absorption 
[22], the  wet  electron  absorption  spectra  is  obtained,  Figure  6.2,  which  is  remarkablely 
similar  to  the  solvated  electron  absorption  spectrum  in  shape,  and  a  red  shift  in 
energy. 
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Figure  6.2 

Wet  Electron  Absorption  Spectra 
wet  electron  (dots)  and  the  solvated  electron  (line) 


Based  on  isosbestic  wavelength  and  wet  electron  absorption  spectrum,  we 
hypothesized  that  the  wet  electron  is  just  the  lowest  excited  state  of  the  solvated 
electron.  Our  experimental  results  show  that  the  wet  electron  formation  time,  1  /ki ,  is 
300  ±  30  fs.  The  lifetime  of  the  wet  electron,  1/k2,  was  found  to  be  540  ±  30  fs.  We 
therefore  see  that  the  rate-limiting  step  in  the  solvation  dynamics  of  the  electror  can 
be  really  an  energy  relaxation  process  like  as  a  radiationless  transition  in  molecules. 
One  anticipates  that  the  internal  conversion  process  can  be  described  by  the  terms 
neglected  in  the  Born-Oppenheimer  approximation,  as  is  the  case  with  molecules.  For 
the  solvated  electron,  it  is  reasonable  to  say  that  the  relevant  motions  might  be  water 
librations.  Recent  calculations  of  Rossky,  Freisner  and  coworkers  provide  theoretical 
support  for  the  excited  state  picture.  [23] 

The  role  of  radiationless  transitions  in  the  solvation  dynamics  of  an  electron 
suggests  that  the  solvated  electron  can  be  viewed  as  a  kind  of  molecule  which 
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undergoes  a  rapid  internal  conversion  to  the  ground  state.  In  the  traditional  picture  of 
molecular  internal  conversion,  Frank-Condon  factors  are  the  dominant  terms  in  the 
golden  rule  expression  for  the  non-radiative  rate. [24]  To  begin  to  address  this 
question,  we  have  examined  the  effects  of  isotope  substitution  on  the  solvation 
dynamics  of  the  electron. 

The  solvated  electron  spectra  in  H20  and  D20  are  the  same  within  a  few 
percent,  thus  it  is  quite  reasonable  to  assume  that  the  wet  electron  absorption  spectra 
are  also  nearly  identical. [22]  Therefore  differences  in  the  observed  solvation 
dynamics  can  be  directly  attributed  to  kinetic  effects.  The  wet  electron  formation  and 
decay  can  be  separated  by  observing  the  dynamics  at  the  isosbestic  wavelength.  The 
dynamics  at  820  nm  are  different  for  H20  and  D20,  1/K1  are  300±30  fs  and  360±30 
fs  respectively.  At  1100  nm  and  625  nm,  clear  differences  are  observed  also.  From 
fitting  the  data  to  equation,  values  of  1/K2,  the  solvated  electron  formation  time,  were 
found  to  be  540±30  fs  and  600±30  fs.  These  results  are  consistent  with  the  overall 
increase  in  the  average  solvation  time  reported  earlier.  The  large  increase  in  the 
solvation  time  upon  isotope  substitution  suggests  that  librational  motions  are  involved 
in  the  rate-limiting  step  in  the  solvation  dynamics. 

6.3  Geminate  Recombination  of  Solvated  Electrons  In  H20  and  D20,  large 
non-exponential  decays  are  seen  in  the  absorption  signal  on  the  picosecond  time 
scale.  All  wavelengths  studied  show  the  same  kinetic  behavior.[25]-  Since  there  is  no 
wavelength  dependence  to  the  data,  it  suggests  strongiy  that  a  single  species  or 
absorber  density  decay  process  is  being  observed.  In  addition  it  is  generally  believed 
that  the  water  cation  is  very  unstable  in  liquid  water  and  within  a  few  tens  of 
femtoseconds  reacts  with  a  neutral  water  to  form  H3O  and  OH. [26]  Thus  we  are 
actually  observing  the  recombination  of  the  solvated  electron  with  the  hydronium 
cation. 

Comparing  the  geminate  recombination  kinetics  of  H20  and  D20,  a  clear 
isotope  effect  is  observed.  We  find  that  the  recombination  is  slower  and  more 
electrons  appear  to  escape  recombination  in  D2O.  To  carefully  verify  the  existence  of 
an  isotope  dependence  on  the  electron  thermalization  distance,  we  have  compared 
our  data  to  the  Onsager  model  of  electron-cation  recombination.  [27]  Results  of  the 
fitting  are  shown  in  Figure  6.3.  Using  an  exponential  distribution  of  solvated  electrons, 
the  thermalization  lengths  of  1 2.4  ±  0.3  A  in  D20  and  1 1 ,5±0.5  A  in  H20  are  obtained 
This  indicates  that  the  isotope  effect  results  in  larger  electron  thermalization  length  in 
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D2O  than  that  in  H2O.  The  physical  basis  of  this  result,  we  believe,  is  the  less  efficient 
energy  transfer  from  the  hot  photoionized  excess  electron  to  the  solvent  modes  in 
D2O  compared  to  H20.[28]  The  reduced  frequencies  mean  that  more  collisions  are 
needed  to  thermalize  the  hot  electron  in  D2O  compared  to  H2O,  thus  results  larger 
thermalization  distance. 


Figure  6.3 

Geminate  Recombination  of  Electrons  in  Water 


The  electron-cation  geminate  recombination  dynamics  studies  were  also 
performed  in  neat  alkane  solvents.  After  photoexcitation  there  can  be  at  least  three 
species  in  the  neat  alkane,  the  electron,  the  alkane  cation  and  the  excited  state  of  the 
alkane.  The  electron  in  an  alkane  absorbs  in  the  infrared,  with  a  peak  around  2000 
nm.  [  e=7000  M^-cnrr1  at  2000  nm  for  n-hexane]  There  is,  however,  a  large 
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absorption  tail  for  the  electron  that  extends  well  into  the  near  IR  and  even  into  the 
visible  at  room  temperature. [29] 

The  absorption  data  for  iso-octane  is  nearly  wavelength  independent.  This 
suggests  that  only  one  single  absorber  or  two  species  with  the  same  kinetics  are  being 
observed.  In  the  IR  it  is  quite  reasonable  to  assume  that  the  electron  will  be  the 
dominant  absorber.  CCI4  quenching  experiments  showed  that  there  is  no  absorption 
of  alkane  cation  at  near  IR.  Linear  alkanes  were  observed  excited  absorption  constant 
backgrounds,  which  can  be  substracted,  at  picosecond  time  scales. 

Qualitatively  there  is  a  clear  correlation  between  the  shape  of  the  molecule,  the 
electron  mobility  and  the  observed  geminate  recombination  kinetics.  Spherical  or 
branched  alkane  molecules  have  comparatively  large  electron  mobilities  and  fast 
geminate  recombination  kinetics.  Linear  alkanes  have  much  smaller  electron 
mobilities  and  exhibit  slow  geminate  recombination  kinetics.  For  example  the 
geminate  recombination  in  iso-octane  is  completed  in  approximately  0.5  ps,  while  the 
geminate  recombination  in  n-octane  is  significant  out  to  50  to  100  ps.  Using  the 
Onsager  model,  fits  were  done  by  varying  mean  thermalization  distance  <r>  of 
exponential  electron  distribution  and  using  the  experimentally  known  values  of 
electron  mobility.  The  above  analysis  yielded  a  mean  electron  thermalization  distance 
of  54±1 1  A  in  iso-octane  and  48±10A  in  normal  octane. 

To  further  examine  the  kinetics  of  geminate  recombination,  the  temperature  of  the  iso¬ 
octane  was  varied.  These  measurements  are  qualitatively  consistent  with  the  diffusion 
picture,  at  higher  temperatures  the  recombination  is  faster.  The  comparison  to  the 
Onsager  model  is  quite  good  and  the  the  thermalization  distances  are  obtained  to  be 
temperature  independent.  This  indicate  that  the  high  frequency  modes  of  solvent 
media  are  mainly  contributed  to  the  kinetic  energy  dissipation  of  excess  electrons. 


3  6 


Relative  Time  ps 


Figure  6.4 
Geminate  Recombination 


Iso-octane 
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Figure  6.5 

Geminate  Recombination  In  n-octane 

6.4  Photodetachment  of  Simple  Aqueous  Anions  Many  simple  anions  in 
solution,  such  as  halides,  have  large  and  broad  absorption  spectra  in  the  UV.  The 
absorption  bands  that  correspond  to  the  removal  of  the  electron  from  the  ion  are  called 
charge  transfer  to  solvent  spectra  or  CTTS  spectra.[301  The  excited  electron  will  only 
experience  a  potential  field  well  due  to  the  orientational  polarization  of  the  water 
molecules,  which  cannot  change  on  the  time  scale  of  optical  excitation. 

The  photodetachment  dynamics  probed  at  625  nm  dramatically  shows  the 
appearance  and  rapid  disappearance  of  a  new  species  in  addition  to  the  wet  and 
solvated  electron. [31]  The  kinetic  behavior  indicates  that  the  species  is  either  the 
electron  or  (|-)\  a  CTTS  excited  state.  The  iodine  atom  is  discounted  because  it  does 
not  absorb  in  the  visible.  It  is  well  established  that  the  excited  state  of  the  the  iodide 
anion  produced  upon  excitation  corresponds  to  the  CTTS  state,  where  the  electron  is 
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trapped  in  the  well  due  to  the  polarization  of  the  water  molecules  around  the  original 
anion. 

Scaning  through  the  visible  and  near  IR  wavelength  range  in  1M  iodide 
aqueous  solution,  we  located  an  isosbestic  point  to  be  at  850  nm.  This  is  very  close  to 
the  820  nm  obtained  for  the  neat  water.  From  a  single  exponential  fit,  the  appearance 
time  of  the  wet  electron  was  found  to  be  500  fs.  This  is  about  50  %  longer  than  the 
time  found  for  neat  water.  However  most  importantly,  the  determination  of  the 
isosbestic  point  shows  that  the  wet  electron  to  solvated  electron  transition  is  a  two- 
state  process  in  nature.  In  difference  from  the  neat  water  case,  the  wet  electron 
precursor  can  be  observed  by  visible  absorption  directly;  the  trapped  electron,  i.e.  the 
electron  in  the  very  short  lived  CTTS  state. 

kl  k2 

©trapped  >©wet  *©sol  (gq.  6.3) 

ki  is  the  formation  rate  of  the  wet  electron;  k2  is  the  formation  rate  of  the  solvated 
electron.  The  data  at  all  wavelengths  were  fitted  to  this  model.  We  obtained  the  wet 
electron  absorption  spectra  in  the  iodide  solution,  by  comparison  with  the  known 
solvated  electron  absorption.  It  is  quite  similar  to  the  wet  electron  photoionized  from  a 
neutral  water  molecule.  The  decay  time  of  the  wet  electron,  l/k2  ,  extracted  from  the 
fit  was  the  same  as  neat  water  540  fs.  This  result  and  the  similarity  of  the  wet  electron 
absorption  spectra  to  that  obtained  for  neat  water  suggest  that  the  largest  deviations 
from  the  neat  water  case  occur  at  early  times.  Clearly  the  situation  is  complex  and 
more  investigation  is  in  progress. 


7.  Picosecond  Dynamics  of  Chemical  Intermediates. 

7.1  Mechanism  of  Singlet  Oxygen  Production  from  the 
Photodissociation  of  Anthracene  Endoperoxides  The  observation  of 
electronic  state-selective  chemistry  provides  valuable  insight  into  the  nature  of 
chemical  reactions.  Aromatic  endoperoxides  are  compounds  known  to  undergo 
wavelength  dependent  photochemistry.  It  is  well  known,  for  example,  that 
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exciting  into  Si  gives  products  derived  from  0-0  rupture,  but  excitation  to 
higher  states  gives  IO2  and  the  anthracene  moiety.  We  extended  our  studies  of 
these  molecules,  characterized  by  an  0-0  peroxide  bridge,  to  better  understand 
the  nature  of  the  photofragmentation  process  and  its  dependence  on  the 
wavelength  of  the  exciting  light.  We  examined  two  anthracene  endoperoxides 
(9,10-DPA-O2  and  9-PA-O2). 


The  mechanism  for  the  state-selective  singlet  oxygen  producing 
photochemistry  of  the  anthracene  endoperoxides  was  determined  using 
picosecond  laser  induced  fluorescence  methods.  For  9.10-DPA-O2  the  ground 
state  aromatic  appears  in  95  ±  10  ps,  whereas  for  9.10-PA-O2  it  appears  in  60  ± 
7  ps.  We  thus  found  that  the  reaction  proceeds  efficiently  from  an  upper  excited 
singlet  state  and  that  it  involves  a  relatively  long  lived  (50  -  100  ps) 
intermediate,  which  was  describable  as  having  significant  biradical  character. 
These  results  are  consistent  with  our  original  interpretation  given  to  the 
photochemistry  of  two  other  anthracene  endoperoxides  (9,10-DMDPA-C>2  and 
1,4-DMDPA-02). 


9,10-DMDPA-O2  1.4-DMDPA-02 

Figure  7.2 

Before  our  work  it  was  widely  assumed  that  photodissociation  of  endoperoxides 
was  concerted,  i.e.,  it  involved  the  simultaneous  rupture  of  both  C-0  bonds 
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leading  directly  to  102  and  the  ground  state  aromatic  moiety.  Consistent  with 
the  observed  kinetics  as  well  as  the  structural  dependence  on  the  rates  is  the 
following  mechanism  for  102  production  given  by  Scheme  1. 

in  summary,  the  observation  of  the  delayed  appearance  of  the 
anthracene  following  photoexcitation  into  the  upper  k,k*  Sn  (n,  >2)  states  as 
well  as  structural  effects  on  the  kinetics  can  be  explained  by  an  intermediate 
that  is  formed  either  by  non-concerted,  sequential  C-0  bond  rupture,  which 
involves  a  1,6-biradical  intermediate,  or  by  the  formation  of  an  excited  state 
anthracene  valence  isomer  that  has  substantial  biradical  character. 


Ph 


Ph 


Figure  7.3 

Biradlcal  mechanism  for  upper  excited  state  n,iz*  state  photofragmentation  of 

aromatic  endoperoxldes. 


7.2  The  Chemical  Intermediate  Dimesityl-carbene  -  Energetics  and 
Intersystem  Crossing  Dynamics 

It  is  well  known  that  derivatives  of  diphenylcarbene  exhibit  significantly 
different  reaction  profiles  toward  various  spin  state  selective  substrates  (such  as 
alcohols  and  olefins).  Based  on  our  work  and  of  others,  we  now  know  that  a  key 
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reason  for  the  variation  in  reaction  efficiencies  for  the  different  carbene 
derivatives  is  related  to  corresponding  differences  in  the  energy  gap  (AEst) 
between  the  lowest  singlet  and  ground  triplet  state  of  the  carbenes.  Changes  in 
AEst  mean  that  singlet/triplet  population  ratios  change  and  thus  will  affect  the 
yields  of  reaction  involving  state  selective  quenchers.  From  our  earlier  work  we 
also  showed  that  solvent  polarity  dramatically  alters  AEst  for  a  given  carbene 
(AEst  decreases  with  increasing  solvent  polarity),  which  affects  not  only  the 
population  of  singlet  and  triplet  states  but  also  the  intersystem  crossing 
dynamics  between  the  two  states. 

A  particularly  atypical  diphenylcarbene  derivative  is  dimesitylcarbene 
(DMC).  For  instance  the  chemisty  of  DMC  appears  to  occur  almost  exclusively 
from  the  ground  triplet 


Dimesitylcarbene  (DMC) 

state,  unlike  that  for  other  diphenylcarbenes.  This  could  be  due  to  an  unusually 
large  AEST  for  DMC,  which  would  make  the  singlet  state  inaccessible  during 
the  carbene  lifetime.  A  key  issue  therefore  is  to  establish  the  energetics  and 
spin  relaxation  dynamics  for  DMC  and  to  relate  it  to  the  other  carbenes. 

Using  picosecond  laser  methods  we  succeeded  in  determining  the  spin- 
relaxation  dynamics  as  well  as  to  establish  a  lower  limit  on  AEst  for  DMC. 
From  our  results  we  were  able  to  answer  the  long  standing  question  about  the 
intersystem  crossing  dynamics  of  dimesitylcarbene  (DMC)  and  how  it  is  related 
to  other  homologous  carbenes.  From  the  solvent  polarity  dependence  we 
deduced  that  intersystem  crossing  in  DMC  proceeds  in  the  large  gap  limit,  in 
which  the  rate  of  singlet-to-triplet  conversion  increases  as  the  energy  gap 
decreases.  This  result  is  consistent  with  the  large  AEst  (>  7-5  kcal/mole)  that 
we  found  for  this  carbene.  The  large  energy  gap  was  explained  by  the  large 
bond  angle  about  the  central  methylene  carbon  atom,  which  is  caused  by  steric 
effects  involving  the  two  mesityl  groups.  This  work  provides  important  new 
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information  about  the  nature  of  the  singlet-triplet  splitting  and  its  relation  to  ihc 
dynamics  associated  with  intersystem  crossing  between  the  lowest  singlet  and 
triplet  states  of  aromatic  carbenes. 

7.3  Photochemistry  in  Restricted  Environments  We  investigated  how 
specific  local  or  heterogeneous  microenvironments  due  to  molecular 
encapsulation  significantly  modify  photoisomerization  reaction  rates.  In  the  case 
of  t-stilbene  complexed  to  various  cyclodextrins  we  directly  observed  how 
frictional  effects  due  to  the  local  (or  "restricted")  environment  of  the  inner 
cyclodextrin  cavity  modifed  the  motion  necessary  for  reaction,  which  thereby 
altered  reaction  rates.  The  effects  caused  by  these  micro-heterogeneous 
environments  could  be  rationalized  by  a  simple  "lock  and  key"  hypothesis, 
which  is  widely  used  to  explain  the  nature  of  the  guest-host  interactions  in 
enzymatic  systems. 

The  trans-cis  photoisomerization  of  stilbene  is  known  to  occur  by  a 
thermally  activated  twisting  around  the  central  double  bond, 


Figure  7.5 

Photoisomerization  of  stilbene 

which  thus  provides  the  principal  mode  of  deactivation  of  the  excited  singlet 
state  of  stilbene.  From  picosecond  time  resolved  methods  we  have  been  able  to 
show  for  the  first  time  that  the  nature  of  the  complex  is  different  for  different 
cyclodextrins,  being  dependent  mainly  on  the  cavity  size  of  the  host  molecule 
(5.6  A  for  a-cyclodextrin  and  up  to  8  A  for  g-cyclodextrin).  This  space  limitation 
affected  the  chemistry  of  a  reactive  guest  molecule  because  it  restricted 
motions  that  were  necessary  for  the  ultrafast  photo-reaction  and  thus  altered  the 
rates  of  reaction.  A  single  complex,  with  simple  reaction  kinetics,  was  found  for 
the  sti'bene/a-cyclodextrin  complex.  This  is  to  be  contrasted  with  the  larger 


AFOSR-88-0014 


"tiaht"  complex 


"loose"  complex 


Figure  7.6 

Schematic  representation  of  the  dynamic  equilibrium  between  "tight"  and 
"loose"  stllbene/cyclodextrln  complexes 


sized  cyclodextrins,  which  yielded  two  types  of  1 :1  complexes  with  the  stilbene 
molecule  (e.g.  weakly  and  tightly  bound  forms).  The  two  distinctly  different 
complexes  formed  for  the  larger  b-  and  g-cyclodextrins,  yielded  vastly  different 
kinetics  for  the  photoisomerization  process.  For  the  "tight"  complex  the  guest 
stilbene  molecule  is  immobilized,  yielding  slow  kinetics  (being  comparable  to 
the  kinetics  observed  for  the  molecule  in  the  solid  state).  In  contrast,  the  "loose” 
complex  undergoes  rapid  photoisomerization,  and  is  explained  by  the  fact  that 
the  local  friction  (viscosity)  experienced  by  the  stilbene  molecule  is  low. 
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List  of  Publications 

"Energetics  of  Adsorption  of  Neutral  and  Charged  Molecules  at  the  Air/Water 
Interface  by  Second  Harmonic  Generation  -  Hydrophobic  and  Solvation  Effects" 
Castro,  A.,  Bhattacharyya,  K.  and  Eisenthal,  K.B. 
accepted  for  publication  in  J.  Chem.  Phys.,  4/91 

"Femtosecond  studies  of  electron  photodetachment  from  an  iodide  in  solution: 
The 

trapped  electron" 

Long,  F.  H.;  Lu,  H.;  Shi,  X.;  Eisenthal,  K.  B. 

Chem.  Phys.  Lett.  169.  165  1990 

"Determination  of  pKa  at  the  air/water  interface  by  second  harmonic  generation” 
Zhao,  X.;  Subrahmanyan,  S.;  Eisenthal.  K.  B. 

Chem.  Phys.  Lett.,  171, 558  (  1990) 

"Femtosecond  studies  of  the  pre-solvated  electron:  An  excited  state  of  the 
solvated  electron?" 

Long,  F.  H.;  Lu,  H.;  Eisenthal,  K.  B. 

Phys.  Rev.  Lett.  M,  1469,  1990 

"Femtosecond  studies  of  electrons  in  liquids" 

Lu,  H.;  Long,  F.  H.;  Eisenthal,  K.  B. 

J.  Opt.  Soc.  Am.  B,7,  1511  {  1990) 

"Studies  of  surface  diffusion  by  second  harmonic  fluctuation  spectroscopy" 
Zhao,  X.;  Goh,  M.  C.;  Subrahmanyan,  S.;  Eisenthal,  K.  B. 

J.  Phys.  Chem.  24,  3370,  1990 

"Effects  of  long  chain  molecules  on  the  population  and  orientaion  of  organic 
molecules 

at  the  air/water  interface" 

Zhao,  X.;  Goh,  M.  C.;  Eisenthal  K.  B. 

J.  Phys.  Chem.  94 , 222  (1990). 

Reply  to  comments  on  "Femtosecond  studies  of  electron-cation  geminate 
recombination  in  water" 

Eisenthal,  Kenneth  B. 

J.  Phys.  Chem.,  93,  7534,  1989 

"Studies  of  molecular  properties  at  the  surface  of  a  liquid  jet  by  second 

harmonic 

generation" 

Castro,  A.;  Ong,  S.;  Eisenthal,  K.  B. 

Chem.  Phys.  Lett.,  163(4,5),  412-416,  1989 

"Femtosecond  studies  of  electron  photodetachment  of  simple  ions  in  liquid 
water: 
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Solvation  and  geminate  recombination  dynamics" 

Long,  F.H.;  Lu,  H.;  Eisenthal,  K.B. 

J.  Chem.  Phys.,  91(7),  4413-4,  1989 

"Effects  of  solvent  polarity  and  structure  on  intersystem  crossing  in 
diphenylcarbenes: 

A  picosecond  laser  study  on  dimesitylcarbene" 

Sitzmann,  E.  V.;  Langan,  J.  G.;  Griller,  D.;  Eisenthal,  K.B. 

Chem.  Phys.  Lett.,  161(4,5),  353-360,  1989 


"Picosecond  kinetics  of  state-selective  singlet  oxygen  producing  photochemistry 
of  aromatic  endoperoxides" 

Sitzmann,  E.  V.;  Langan,  J.  G.;  Hrovat,  D.  A.;  Eisenthal,  K.  B. 

Chem.  Phys.  Lett.,  162(1,2),  157-162,  1989 

"Femtosecond  studies  of  electron-cation  dynamics  in  neat  water:  The  effects  of 
isotope  substitution" 

Long,  Frederick  H.;  Lu,  Hong;  Eisenthal,  Kenneth  B. 

Chem.  Phys.  Lett.,  160(4),  464-8,  1989 

"The  role  of  translational  friction  in  isomerization  reactions" 

Bowman,  Robert  M.;  Eisenthal,  Kennth  B. 

Chem.  Phys.  Lett.,  155(1),  99-101,  1989 

"Picosecond  laser  studies  on  photochemical  reactions  in  restricted 
environments: 

The  photoisomerization  of  t  -stilbene  complexed  to  cyclodextrins" 

Duveneck,  G.L.;  Sitzmann,  E.V.;  Eisenthal,  K.B.;  Turro,  N.J. 

J.  Phys.  Chem.,  93(20),  7166-7170,  1989 

"Dynamics  of  intermolecular  electronic  energy  transfer  at  an  air/liquid  interface" 
Sitzmann,  Eugene  V.;  Eisenthal,  Kenneth  B. 

J.  Chem.  Phys.,  90(5),  2831-2,  1989 

"The  energetics  of  orientation  at  the  liquid-vapor  interface  of  water" 

Goh,  M.  Cynthia;  Eisenthal,  Kenneth  B. 
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Invited  Talks  1989-1991 

American  Chemical  Society  Symposium  on  "Energy  Transfer  and 
Relaxation,"  25-30  August,  1991.  Invited  Speaker. 

EUCMOS  XX,  25-30  August  1991,  Zagreb,  Yugoslavia.  Invited  Speaker. 

International  Symposium  on  Ultrafast  Processes  in  Spectroscopy, 
October  1991.  Bayreuth,  Germany.  Invited  Speaker. 

Symposium  on  "Dynamical  Processes  in  Condensed  Molecular 

Systems,  29-30  April  1991.  Neve-Ilan,  Israel.  Invited  Speaker. 

International  Conference  on  Coherent  and  Nonlinear  Optics, 

24-27  September  1991,  Leningrad,  USSR.  Invited  Speaker. 

International  Workshop  on  Laser  Physics,  Dalian  Institute  of 
Chemical  Physics,  19  May  -  1  June  1991.  Dalian,  China. 

Invited  Speaker. 

Center  for  Analysis  of  Structures  and  Interfaces  Colloquium  "Recent 
Advances  in  the  Uses  of  Light  in  Physics,  Chemistry  and 
Medicine,  City  College,  CUNY.  19-21  June,  1991.  Invited 
Speaker. 

American  Physical  Society  Symposium  on  Interfaces,  18-22  March 
1991.  Cincinnati,  Ohio.  Invited  Speaker. 

XV-th  International  Conference  on  Photochemistry,  28  July  - 
2  August  1991,  Paris,  France.  Invited  Speaker. 


American  Chemical  Society,  Symposium  on  Classical  and  Quantal 
Simulations  for  Reactive  and  Solvation  Dynamics  and  Their 
Critical  Experimental  Tests,  23-27  April  1990,  Boston, 
Massachusetts.  Invited  speaker. 

Argonne  National  Laboratory,  Symposium  on  The  Solvated 

Electron-Past,  Present  and  Future,  16-17  July  1990,  Argonne, 
Illinois.  Invited  speaker. 
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Gordon  Research  Conference  on  the  Chemistry  and  Physics  of  Water 
and  Aqueous  Solutions,  30  July-3  August  1990,  Plymouth,  New 
Hampshire.  Invited  speaker. 

Gordon  Research  Conference  on  Radiation  Chemistry,  9-13  July  1990, 
Salve  Regina,  Rhode  Island.  Invited  speaker. 

International  Centre  for  Theoretical  Physics,  Conference  on  Lasers  in 
Chemistry,  11-15  June  1990,  Miramare-Trieste,  Italy.  Invited 
speaker. 

Naval  Research  Laboratory,  9  March  1990,  Washington,  D.  C.  Invited 
speaker. 

Optical  Society  of  America,  Seventh  Topical  Meeting  on  Ultrafast 
Phenomena,  14-17  May  1990,  Monterey,  California.  Invited 
speaker. 

Research  Foundation  for  Opto-science  and  Technology,  International 
Seminar  on  Photochemical  Dynamics,  25-27  March  1990, 
Hamamatsu,  Japan.  Invited  speaker. 

State  University  of  New  York,  Buffalo,  21  March  1990.  Invited 
speaker. 

University  of  Michigan,  Department  of  Atmospheric,  Oceanic,  and 
Space  Sciences,  Nineteenth  Informal  Conference  on 
Photochemistry,  25-29  June  1990,  Ann  Arbor,  Michigan. 

Invited  speaker. 

U.S.A.-U.S.S.R.  Joint  Workshop  on  Linear  and  Nonlinear  Laser 
Interactions  and  Molecular  Dynamics,  18-22  June  1990, 

Moscow,  Leningrad,  U.S.S.R.  Invited  speaker. 

Rochester  American  Chemical  Society  Section’s  Fall  Symposium, 
University  of  Rochester,  October  8,  1990.  Invited  Speaker. 

NATO  Advanced  Research  Workshop,  Tutzing,  Bavaria, 

10-13  September,  1990.  Discussion  Leader. 

Naval  Research  Laboratory,  Washington,  D.C.,  9  March  1990. 

Invited  Speaker. 
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University  of  Colorado,  Boulder,  16  March  1990.  Invited  Speaker. 

State  University  of  New  York  at  Buffalo,  21  March  1990.  Invited 
Speaker. 

Union  Carbide  Corporation,  Tarrytown,  N.Y.,  23  May,  1990.  Invited 
Speaker. 

University  of  Wisconsin,  Madison.  September,  1990.  Invited 
Speaker. 

AT&T  Bell  Laboratories,  Murray  Hill,  N.J.,  23  October,  1990.  Invited 
Speaker. 

American  Chemical  Society,  Division  of  Physical  Chemistry, 

Femtosecond  Chemistry  Symposium,  10-15  September  1989, 
St.  Louis,  Missouri.  Invited  speaker. 

American  Chemical  Society,  Division  of  Physical  Chemistry, 
Symposium  on  Chemistry  at  Solution  Interfaces, 

9-14  April  1989,  Dallas,  Texas.  Invited  speaker. 

American  Physical  Society  Meeting,  Symposium  of  the  Division  of 
Chemical  Physics:  Reaction  Dynamics  in  Liquids  on  Ultrafast 
Time  Scales,  20-24  March  1989,  St.  Louis,  Missouri. 

Invited  speaker. 

California  Institute  of  Technology,  18  December  1989.  Invited 
speaker. 

Gordon  Research  Conference  on  Liquids,  13-18  August  1989, 
Plymouth,  New  Hampshire. 

Optical  Society  of  America,  Quantum  Electronics  and  Laser  Science 
Conference,  24-28  April  1989,  Baltimore,  Maryland.  Invited 
speaker. 

The  Pennsylvania  State  University,  19  September  1989.  Invited 
speaker. 

Polytechnic  University,  27  September  1989.  Invited  speaker. 
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Princeton  University,  Fourth  International  Conference  on  Time- 

resolved  Vibrational  Spectroscopy,  11-16  June  1989,  Princeton, 
New  Jersey.  Invited  speaker. 

Stanford  University,  5  December  1989.  Invited  speaker. 

University  of  California,  Berkeley,  21  November  1989.  Incited 
speaker. 

University  of  California,  Irvine,  11  December  1989.  Invited  speaker. 
University  of  California,  San  Diego,  14  March  1989.  Invited  Speaker. 
University  of  Rochester,  4  October  1989.  Invited  speaker. 


